Benchtop Models of Patient-Specific Intraventricular Flow During Heart Failure and LVAD Support.

The characterization of intraventricular flow is critical to evaluate the efficiency of fluid transport and potential thromboembolic risk but challenging to measure directly in advanced heart failure (HF) patients with left ventricular assist device (LVAD) support. The study aims to validate an in-house mock loop (ML) by simulating specific conditions of HF patients with normal and prosthetic mitral valves (MV) and LVAD patients with small and dilated left ventricle volumes, then comparing the flow-related indices result of vortex parameters, residence time (RT), and shear-activation potential (SAP). Patient-specific inputs for the ML studies included heart rate, end-diastolic and end-systolic volumes, ejection fraction, aortic pressure, E/A ratio, and LVAD speed. The ML effectively replicated vortex development and circulation patterns, as well as RT, particularly for HF patient cases. The LVAD velocity fields reflected altered flow paths, in which all or most incoming blood formed a dominant stream directing flow straight from the mitral valve to the apex. RT estimation of patient and ML compared well for all conditions, but SAP was substantially higher in the LVAD cases of the ML. The benchtop system generated comparable and reproducible hemodynamics and fluid dynamics for patient-specific conditions, validating its reliability and clinical relevance. This study demonstrated that ML is a suitable platform to investigate the fluid dynamics of HF and LVAD patients and can be utilized to investigate heart-implant interactions.

Non-Newtonian blood rheology impacts left atrial stasis in patient-specific simulations.

The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects are subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered when calculating patient-specific blood stasis indices by computational fluid dynamics.

Demonstration of Patient-Specific Simulations to Assess Left Atrial Appendage Thrombogenesis Risk.

Atrial fibrillation (AF) alters left atrial (LA) hemodynamics, which can lead to thrombosis in the left atrial appendage (LAA), systemic embolism and stroke. A personalized risk-stratification of AF patients for stroke would permit improved balancing of preventive anticoagulation therapies against bleeding risk. We investigated how LA anatomy and function impact LA and LAA hemodynamics, and explored whether patient-specific analysis by computational fluid dynamics (CFD) can predict the risk of LAA thrombosis. We analyzed 4D-CT acquisitions of LA wall motion with an in-house immersed-boundary CFD solver. We considered six patients with diverse atrial function, three with either a LAA thrombus (removed digitally before running the simulations) or a history of transient ischemic attacks (LAAT/TIA-pos), and three without a LAA thrombus or TIA (LAAT/TIA-neg). We found that blood inside the left atrial appendage of LAAT/TIA-pos patients had marked alterations in residence time and kinetic energy when compared with LAAT/TIA-neg patients. In addition, we showed how the LA conduit, reservoir and booster functions distinctly affect LA and LAA hemodynamics. Finally, fixed-wall and moving-wall simulations produced different LA hemodynamics and residence time predictions for each patient. Consequently, fixed-wall simulations risk-stratified our small cohort for LAA thrombosis worse than moving-wall simulations, particularly patients with intermediate LAA residence time. Overall, these results suggest that both wall kinetics and LAA morphology contribute to LAA blood stasis and thrombosis.

Intraventricular Flow Patterns in Patients Treated with Left Ventricular Assist Devices.

The success of left ventricular assist device (LVAD) therapy is hampered by complications such as thrombosis and bleeding. Understanding blood flow interactions between the heart and the LVAD might help optimize treatment and decrease complication rates. We hypothesized that LVADs modify shear stresses and blood transit in the left ventricle (LV) by changing flow patterns and that these changes can be characterized using 2D echo color Doppler velocimetry (echo-CDV). We used echo-CDV and custom postprocessing methods to map blood flow inside the LV in patients with ongoing LVAD support (Heartmate II, N = 7). We compared it to healthy controls (N = 20) and patients with dilated cardiomyopathy (DCM, N = 20). We also analyzed intraventricular flow changes during LVAD ramp tests (baseline ± 400 rpm). LVAD support reversed the increase in blood stasis associated with DCM, but it did not reduce intraventricular shear exposure. Within the narrow range studied, the ventricular flow was mostly insensitive to changes in pump speed. Patients with significant aortic insufficiency showed abnormalities in blood stasis and shear indices. Overall, this study suggests that noninvasive flow imaging could potentially be used in combination with standard clinical methods for adjusting LVAD settings to optimize flow transport and minimize stasis on an individual basis.

Blood Stasis Imaging Predicts Cerebral Microembolism during Acute Myocardial Infarction.

BACKGROUND: Cardioembolic stroke is a major source of mortality and disability worldwide. The authors hypothesized that quantitative characterization of intracardiac blood stasis may be useful to determine cardioembolic risk in order to personalize anticoagulation therapy. The aim of this study was to assess the relationship between image-based metrics of blood stasis in the left ventricle and brain microembolism, a surrogate marker of cardiac embolism, in a controlled animal experimental model of acute myocardial infarction (AMI). METHODS: Intraventricular blood stasis maps were derived from conventional color Doppler echocardiography in 10 pigs during anterior AMI induced by sequential ligation of the mid and proximal left anterior descending coronary artery (AMI-1 and AMI-2 phases). From these maps, indices of global and local blood stasis were calculated, such as the average residence time and the size and ratio of contact with the endocardium of blood regions with long residence times. The incidence of brain microemboli (high-intensity transient signals [HITS]) was monitored using carotid Doppler ultrasound. RESULTS: HITS were detected in 0%, 50%, and 90% of the animals at baseline and during AMI-1 and AMI-2 phases, respectively. The average residence time of blood in the left ventricle increased in parallel. The residence time performed well to predict microemboli (C-index = 0.89, 95% CI, 0.75-1.00) and closely correlated with the number of HITS (R = 0.87, P < .001). Multivariate and mediation analyses demonstrated that the number of HITS during AMI phases was best explained by stasis. Among conventional echocardiographic variables, only apical wall motion score weakly correlated with the number of HITS (R = 0.3, P = .04). Mural thrombosis in the left ventricle was ruled out in all animals. CONCLUSIONS: The degree of stasis of blood in the left ventricle caused by AMI is closely related to the incidence of brain microembolism. Therefore, stasis imaging is a promising tool for a patient-specific assessment of cardioembolic risk.

Regional dynamics of fractal dimension of the left ventricular endocardium from cine computed tomography images.

We present a method to leverage the high fidelity of computed tomography (CT) to quantify regional left ventricular function using topography variation of the endocardium as a surrogate measure of strain. 4DCT images of 10 normal and 10 abnormal subjects, acquired with standard clinical protocols, are used. The topography of the endocardium is characterized by its regional values of fractal dimension ( F D ), computed using a box-counting algorithm developed in-house. The average F D in each of the 16 American Heart Association segments is calculated for each subject as a function of time over the cardiac cycle. The normal subjects show a peak systolic percentage change in F D of 5.9 % ± 2 % in all free-wall segments, whereas the abnormal cohort experiences a change of 2 % ± 1.2 % ( p < 0.00001 ). Septal segments, being smooth, do not undergo large changes in F D . Additionally, a principal component analysis is performed on the temporal profiles of F D to highlight the possibility for unsupervised classification of normal and abnormal function. The method developed is free from manual contouring and does not require any feature tracking or registration algorithms. The F D values in the free-wall segments correlated well with radial strain and with endocardial regional shortening measurements.

Age-Dependence of Flow Homeostasis in the Left Ventricle.

Background: Intracardiac flow homeostasis requires avoiding blood stasis and platelet activation during its transit through the cardiac chambers. However, the foundations of intraventricular blood washout and its exposure to shear stresses have been poorly addressed. We aimed to characterize and quantify these features in a wide population of healthy subjects and assess the relationships of these indices with age. Methods: We used color-Doppler echocardiography and custom post-processing methods to study 149 healthy volunteers from 26 days to 80 years old. From the intraventricular flow-velocity fields we obtained personalized maps of (1) the residence time of blood in the LV, and (2) the shear index, a metric accounting for the strongest occurrence of shear stresses inside the chamber. From these maps we derived quantitative indices of the overall intraventricular blood washout and shear exposure. We addressed the age-dependence of these indices and analyzed their relationship with age-related changes in filling-flow. Results: The entire intraventricular blood pool was replaced before 8 cycles. Average residence time of blood inside the LV was <3 cycles in all subjects and followed an inverse U-shape relationship with age, increasing from median (IQR) of 1.0 (0.7 to 1.2) cycles in the 1st year of life to 1.8 (1.4-2.2) cycles in young adults (17-30 years old), becoming shorter again thereafter. Shear index showed no relation with age and was bounded around 20 dyn·s/cm2. Regions with the longest residence time and highest shear index were identified near the apex. Differences in the degree of apical penetration of the filling waves and the duration of the late-filling phase explained the age-dependence of residence time ( R adj 2 = 0.48, p < 0.001). Conclusions: In average, blood spends 1 to 3 beats inside the LV with very low shear stress rates. The apical region is the most prone to blood stasis, particularly in mid-aged adults. The washout of blood in the normal LV is age-dependent due to physiological changes in the degree of apical penetration of the filling waves.

Mitral Valve Prosthesis Design Affects Hemodynamic Stasis and Shear In The Dilated Left Ventricle.

Dilated cardiomyopathy produces abnormal left ventricular (LV) blood flow patterns that are linked with thromboembolism (TE). We hypothesized that implantation of mechanical heart valves non-trivially influences TE risk in these patients, exacerbating abnormal LV flow dynamics. The goal of this study was to assess how mitral valve design impacts flow and hemodynamic factors associated with TE. The mid-plane velocity field of a silicone dilated LV model was measured in a mock cardiovascular loop for three different mitral prostheses, two with multiple orientations, and used to characterize LV vortex properties through the cardiac cycle. Blood residence time and a platelet shear activation potential index (SAP) based on the cumulative exposure to shear were also computed. The porcine bioprosthesis (BP) and the bileaflet valve in the anti-anatomical (BL-AA) position produced the most natural flow patterns. The bileaflet valves experienced large shear in the valve hinges and recirculating shear-activated flow, especially in the anatomical (BL-A) and 45-degree (BL-45) positions, thus exhibited high SAP. The tilting disk valve in the septal orientation (TD-S) produced a complete reversal of flow and vortex properties, impairing LV washout and retaining shear-activated fluid, leading to the highest residence time and SAP. In contrast, the tilting disk valve in the free-wall position (TD-F) exhibited mid-range values for residence time and SAP. Hence, the thrombogenic potential of different MHV models and configurations can be collectively ranked from lowest to highest as: BP, BL-AA, TD-F, BL-A, BL-45, and TD-S. These findings provide new insight about the effect of fluid dynamics on LV TE risk, and suggest that the bioprosthesis valve in the mitral position minimizes this risk by producing more physiological flow patterns in patients with dilated cardiomyopathy.

Clinical assessment of intraventricular blood transport in patients undergoing cardiac resynchronization therapy.

In the healthy heart, left ventricular (LV) filling generates different flow patterns which have been proposed to optimize blood transport by coupling diastole and systole. This work presents a novel image-based method to assess how different flow patterns influence LV blood transport in patients undergoing cardiac resynchronization therapy (CRT). Our approach is based on solving the advection equation for a passive scalar field from time-resolved blood velocity fields. Imposing time-varying inflow boundary conditions for the scalar field provides a straightforward method to distinctly track the transport of blood entering the LV in the different filling waves of a given cardiac cycle, as well as the transport barriers which couple filling and ejection. We applied this method to analyze flow transport in a group of patients with implanted CRT devices and a group of healthy volunteers. Velocity fields were obtained using echocardiographic color Doppler velocimetry, which provides two-dimensional time-resolved flow maps in the apical long axis three-chamber view of the LV. In the patients under CRT, the device programming was varied to analyze flow transport under different values of the atrioventricular conduction delay, and to model tachycardia (100 bpm). Using this method, we show how CRT influences the transit of blood inside the left ventricle, contributes to conserving kinetic energy, and favors the generation of hemodynamic forces that accelerate blood in the direction of the LV outflow tract. These novel aspects of ventricular function are clinically accessible by quantitative analysis of color-Doppler echocardiograms.

A clinical method for mapping and quantifying blood stasis in the left ventricle.

In patients at risk of intraventrcular thrombosis, the benefits of chronic anticoagulation therapy need to be balanced with the pro-hemorrhagic effects of therapy. Blood stasis in the cardiac chambers is a recognized risk factor for intracardiac thrombosis and potential cardiogenic embolic events. In this work, we present a novel flow image-based method to assess the location and extent of intraventricular stasis regions inside the left ventricle (LV) by digital processing flow-velocity images obtained either by phase-contrast magnetic resonance (PCMR) or 2D color-Doppler velocimetry (echo-CDV). This approach is based on quantifying the distribution of the blood Residence Time (TR) from time-resolved blood velocity fields in the LV. We tested the new method in illustrative examples of normal hearts, patients with dilated cardiomyopathy and one patient before and after the implantation of a left ventricular assist device (LVAD). The method allowed us to assess in-vivo the location and extent of the stasis regions in the LV. Original metrics were developed to integrate flow properties into simple scalars suitable for a robust and personalized assessment of the risk of thrombosis. From a clinical perspective, this work introduces the new paradigm that quantitative flow dynamics can provide the basis to obtain subclinical markers of intraventricular thrombosis risk. The early prediction of LV blood stasis may result in decrease strokes by appropriate use of anticoagulant therapy for the purpose of primary and secondary prevention. It may also have a significant impact on LVAD device design and operation set-up.